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Abstract. The Einstein/ Abelian-Yang-Mills Equations reduce in the 
^ • stationary and axially symmetric case to a harmonic map with pre- 

D I scribed singularities (p: R^\E ^ IHlJ;^^ into the (fc + l)-dimensional com- 

^0 ' plex hyperbohc space. In this paper, we prove the existence and unique- 

ness of harmonic maps with prescribed singularities ip : R" \ E H, 
where E is an unbounded smooth closed submanifold of R" of codi- 
mension at least 2, and H is a real, complex, or quaternionic hyper- 
bolic space. As a corollary, we prove the existence of solutions to the 
^ , reduced stationary and axially symmetric Einstein/ Abelian-Yang-Mills 

CO ' Equations. 

o 
o 

o\ 

. It is well known that the Einstein/Maxwell Equations reduce in the sta- 

Q^ I tionary and axially symmetric case to a harmonic map with prescribed 

■ singularities \ E — > into the complex hyperbolic plane. Here 

W), S = A \ U^i closed set obtained by removing N open intervals Ij 

bJO, from the z-axis A C M^. More generally, the Einstein/ Abelian-Yang-Mills 

' Equations with Abelian gauge group T'^ reduce in the same setting to a har- 

monic map with prescribed singularities \ S ^ I^c"''^' '^^^^re EI^"''^ is 

^ ! the {k + l)-dimensional complex hyperbolic space. This model has k 'elec- 

^ ' tromagnetic' fields which interact weakly via the gravitational field. The 

number corresponds to the number of black holes. 

In [We4], we proved the existence and uniqueness of harmonic maps with 



1. Introduction 



prescribed singularities (p: M"" \ S — > H, where EI is any hyperbolic space 
(real, complex, or quaternionic), and S is any smooth closed submanifold of 
M" of codimension at least two, provided a harmonic map (p : R" \S M was 
given, prescribing the behavior of (p at infinity in M"\S. Here S denotes the 
union of the unbounded components of S. The singular behavior of (p near 
each component of E was prescribed by a 'Ej -singular map into jj, that is 
a harmonic map <pj : M" \ ^ IHI whose image lies within a geodesic of H, 
and such that ^Pj{x) 7j(-|-oo) as x — > T,j. Such maps are easily constructed 
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using harmonic functions. The map ip was required to be asymptotic to 
near Sj, and asymptotic to at infinity. In particular, it followed that the 
map ip was to be asymptotic near each unbounded component of S to 
'Pj . This map (p was to be considered part of the data, and was to be found 
before the main theorem in [We4] could be applied. 

For the harmonic map problem associated with the Einstein/Maxwell 
Equations, the map (p could be obtained from an explicit solution, the Kerr- 
Newman solution Q, corresponding to = 1 black holes. However, the 
existence of the harmonic map ip in the more general situation of the main 
theorem in [We4] is not known. In this paper, we replace the hypothesis on 
the existence of the harmonic map (p with a considerably weaker hypothesis: 
the existence of a map ip which is only approximately harmonic at infinity, 
see the Main Theorem in Section ^ Such maps are much easier to construct 
than harmonic maps. As a corollary, we obtain the existence of solutions 
to the reduced stationary and axially symmetric Einstein/ Abelian- Yang- 
Mills Equations corresponding to ^ 1 black holes, i.e., the existence of 
a harmonic map with prescribed sing ularities ip:M.^\T.^ Hc+\ where S 
is as in the first paragraph of this introduction. Of course, when k = 1, we 
recover our previous result from [We4]. Since we are interested in unbounded 
singular sets S C M" of codimension at least 2, we will assume throughout 
that n ^ 3. 

It is important to point out that, as before in | We2 , We4 |, unless = 
1, the solutions of the Einstein/ Abelian- Yang-Mills Equations constructed 
in this way, are expected to have a conical singularity on some bounded 
component of the axis. This conical singularity can be interpreted as the 
gravitational force between the black holes, see [ BW , We2 |. However, when 

= 1, these spacetimes should be regular, and asymptotically flat. To 
establish this fact, a regularity result as in |Wel, LTi2 | is needed. We will 
come back to this problem in a future paper. 

In the next section, we describe briefly the setup of the problem. For 
more details, we refer the reader to | [We3 , We4]. In section ^, we state our 
main existence result and give its proof. 



2. Set-Up and Preliminaries 

Let n ^ 3, let C M" be a smooth domain, and let X be a Riemannian 
manifold of dimension m. A harmonic map is a map (/?: ^ X which is for 
each 0' CC O a critical point of the energy functional: 

En'{p) = / \dpf , 
Jn' 

where \dp\'^ = X]j=i(^jV'' is the energy density of p, and (•,-)x rep- 
resents the metric on X. If p is harmonic, it satisfies an elliptic system of 
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nonlinear partial differential equations which can be written in local coordi- 
nates on X as: 

m n 
b,c=l j=l 

where T^^ are the Christoffel symbols of X. Note that for any map the 
quantities t"" are the components of a cross section T{ip) of the pull-back 
bundle (p^^Tli. This vector field along ip is called the tension of ip, and its 
magnitude, defined by {rf = (r, r)x, is a real function on Q. 

Here we will be concerned only with harmonic maps into a classical glob- 
ally symmetric space of noncompact type and rank one, i.e., a real, complex, 
or quaternionic hyperbolic space. We denote such a space by EI and let m 
be its real dimension. The model we use for EI is = M x M'""-'^ where the 
first coordinate is a Busemann function on IH, see [EO, WeS]. The metric 
of EI then takes the form: 

(1) ds'^ = du^ + Qp{dv), 

where for each p G EI, Qp is a quadratic form on R™"^^. 

Example 1. For example, when EI is the complex hyperbolic space EI^^^, we 
have m = 2k -\- 2, and we write p = {u,v,X:'^): where X^'^ ^ l^'^- The 
quadratic form Qp is then given by: 

Qp{dv) = e^''{dv + x-diP-iP- dxf + e'" {dx ■ dx + dtjj ■ dijj) , 

where the dot represents the Euclidean inner product on M^, see [ |We3 , 
Lemma 6]. We also give here for future reference the components of the 
tension t((/?) of a map ip = (n, v, -0) into EI^"''^: 

r" = Au - 2e^" | + % • V-0 - -0 • Vx? 

- e^" {Vx • Vx + V-0 • ViA) 

= 4 Vn - 26^'' (% • V?/' + "0 • Vx) Vu 

- 2^" (X • Vx + -0 • ViA) (Vi; + X • VV - 0Vx) 

r'^ = Vx + 2 Vx Vn - 2e^"V0 (V^; + x • Vi/; - • Vx) 

r-^ = V-i/j + 2 V0 Vu + 2e2"Vx (Vw % • Vi/j - -0 • Vx) ■ 

Harmonic maps have been studied extensively. When the target space is 
nonpositively curved, as is the case here, finite energy harmonic maps are 
well understood, see |ES, H, SUl , SU2 |. In particular, if ^p is any finite energy 
harmonic map then (p is smooth. In this paper, we deal with infinite energy 
harmonic maps into IH with prescribed singularities along a submanifold S 
of codimension at least two. We now give a few definitions taken from [ We3 , 
Wei . 

Definition 1. Let Vt C M" be a smooth domain, and let S be a smooth 
closed submanifold of M" of codimension at least two. Let 7 be a geodesic 
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of H. We say that a harmonic map ip: r2\S^EIisa H-singular map into 
7 if: 

(i) (^(1^\S)C7(M); 

(ii) ip{x) — > 7(+oo) as x — > S; 

(iii) There is a constant 5 > such that \d(p{x)\^ ^ 5dist(x,S)^^ in some 
neighborhood of S. 

Note that if u is a harmonic function on 17 \ S then = 7 o n is a harmonic 
map and \dip\ = |Vn| . Thus, the problem of finding S-singular maps into 
7 is simply a problem in potential theory. 

Example 2. Let C and let S C be a smooth closed submanifold of 
codimension at least two. Let uq be the potential of a charge distribution 
which is positive and bounded away from zero on S, and let vq € M™~^ be 
constant. Then, using the coordinate system for EI described above in which 
the metric takes the form (|l]), we see that ip = [uq^vq) is a S-singular map 
into the geodesic ^{t) = (t, vq) of H. Such maps will be used to prescribe 
the singular behavior near S. 

The next definition generalizes the concept of asymptotic geodesies, see 
for example |[EQ(| . 

Definition 2. Let ip,ip' : Q\Ti —>■ M he maps, and let S' C S. We say that 
ip and ip' are asymptotic near S' if there is a neighborhood SI' of S' such 
that dist(^, (^') E L°°(17' \ S'). 

To ensure uniqueness, we will use the following notion for being asymp- 
totic at infinity. 

Definition 3. Let (/?, (p' : M" \ S ^ H be maps. We say that (p and cp' are 

asymptotic at infinity, if dist((/3, ip') — > as x ^ 00 in R" \ S. 

Consider now the Einstein / Abelian- Yang-Mills Equations with an Abelian 
gauge group T^: 

RiCg -^Rgg = 2Tf, 

(3) ^ = dA, 

d*F = 0, 

Tf{X,Y) =ilg{ixF -iyF + ix*F -iy *F). 

Here g \s a, Lorentzian metric on a 4-dimensional manifold M, Ric^ is the 
Ricci curvature of g, Rg its scalar curvature, A is an M'^-valued one-form 
on M, d is the exterior derivative, ix inner multiplication by the vector X, 
iig is the trace with respect to the metric g, and the dot represents a Eu- 
clidean inner product on R^. We impose the condition that the solution is 
asymptotically flat, stationary and axially symmetric, globally hyperbolic, 
and that its domain of outer communications is simply connected. Further- 
more, we will assume that (M, g) has an event horizon with N ^ 1 connected 



components, each of which is nondegenerate, see [We4] for the definitions 
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Following step by step the reduction in | We4 , Section 2], one sees that 



the equations reduce to a harmonic map with prescribed singularities into 
IH^"^"^. For this, one only needs to note that the one-forms a = i^F, and 
(3 = i^*F are now M'^'-valued. Here ^ is the generator of the axial symmetry 
normalized so that its orbits have parameter length 27r. Thus the potentials 
X, and ■0, given by dx = ct, and dij^ = (3 are also M'^- valued. From here on, 
all the calculations go through, with all products between M'^-valued objects 
being taken in the sense of the Euclidean inner product on M'^. In particular, 
the twist of ^ is given by a; = *{d£^ A where * is the Hodge star operator, 
and the twist potential v is defined hy 2 dv = lo — 2[x ■ dxjj — xj) ■ dx)- Define 
= - log 1^1 on M' = {g{C,0 > 0}> and let Q be the quotient of M' by its 
group M X SO{2) of isometrics. Then Q x S0{2) is diffeomorphic to \ S, 
where S = A \ Uj=i^j, Ij 

are open intervals of the z-axis A C M^, and 
is the number of connected components of the event horizon. Let r be 
another generator of the group of isometrics, normalized so that |rp — > —1 
at spacelike infinity in (M, g), and let p be the area element of the orbit, i.e., 
= — A r| . Then Ap = 0, hence /) is a harmonic coordinate on Q, and 
we can pick a conjugate harmonic coordinate z, so that the metric of Q is 
conformal to dp'^ + dz'^. We now put on Q x SO (2) = \ S the flat metric 
dp'^ + dz'^ + p^dcjP'. Let if = (u, ■0), then if is an axially symmetric 
harmonic map into ]H^+^ defined on \ s. Thus, we obtain the following 



result analogous to Theorem 1 in [We4]. 



Theorem 1. Let (M, g, A) he a solution of the Einstein/ Ahelian- Yang- Mills 
Equations with Ahelian gauge group , which is stationary and axially sym- 
metric. Assume that M is simply connected, that {M,g) is the domain of 
outer communications of an asymptotically flat, globally hyperbolic space- 
time, and that the event horizon in (M, g) has N ^ 1 connected components 
each of which is nondegenerate. Define S C and ip: \ S ^c^^ 
above. Then ip is an axially symmetric harmonic map. 

Furthermore the axis regularity conditions imply that the map f is as- 
ymptotic near each T,j to a Sj-singular harmonic map ipj into a geodesic 
7j of EI^^^. To prescribe the behavior at infinity, we now construct a map 
Ip: \ E ^ ^c^^ which is nearly harmonic at infinity. We note that essen- 
tially the same construction will give maps from R" \ S into quaternionic 
hyperbolic space, provided S is of codimension 2, and each component of S 
is contained within some cone centered at the origin. 

Example 3. Let S C R^ be as above. For each 1 ^ j ^ N -\- 1, \et (pj : \ 
S m^"*"^ be a given Sj-singular map into a geodesic jj of H^'*'^. We 
assume, without loss of generality, that all the geodesies 7j have a common 
initial point 7j( — oo) on aEI^+\ so that the maps (pj can all be written as 
in Example |2[ i.e., ipj = {uQ,Vj,Xj-,'^j) where Vj G M, and Xj^'4'j ^ are 
constants. We also assume that uq is the potential of a uniform unit charge 
on S. We can then normalize uq so that uq — »■ — log p as r — > oo in R^ \ S. 
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We seek a map : \ S ^c^'^ which coincides in a neighborhood of each 
Tij with ipj , and which is nearly harmonic at infinity. More precisely, we will 
obtain a map for which there is a constant c > such that 



(4) |r((^)| ^c(l+r2)- 



where r is the distance to the origin in M.^. To achieve this, we simply set 
u = uq, and we pick functions v, x, which depend only on the polar angle 
outside a sufficiently large ball, and which take on the appropriate constant 
values Vj , Xj , V'j ) near each . It is now a simple calculation using (§) to 
check that the map (p = {u,v,X-,'4^) satisfies @. 

We can now pose the following Reduced Problem for the stationary and 
axially symmetric Einstein/ Abelian- Yang-Mills Equations. 

Reduced Problem . Let E = UjS^ = -'^ \ Uf=i where A is the z- 
axis in M^, Ij are open intervals, and Sj are the connected components 
of S. Let Uq be the potential of a uniform unit charge distribution on S. 
For each 1 ^ j ^ N + 1, let vj € M and ifyj G M.^ be constants, and let 
ipj : \ T, ^c^^ be the S-singular map given, as in Example ^, by 
= iu-OjVjjO, il^j)- Construct the map ^p: \ S — >• EI^"^^ as in Example ^. 
The Reduced Problem is then: to find an axially symmetric harmonic map 
\ S — > H^"''^ which is asymptotic to ^ near each Sj, and at infinity. 

Let S' = Si U Siv+i, and denote by p' : \ S' ^ the harmonic map 



obtained from the Kerr-Newman solution |]C|, which was used in [We4] to 
prescribe the behavior of (/? at infinity. Note that and p' are asymptotic 
at infinity. Thus, using ip' is equivalent in this case to using ip. 

Li the next section, we prove a theorem which implies that the Reduced 
Problem has a unique solution for each value of the (A^-|-1)(A;-|-1) parameters. 
We note that one of the constants Vj and one of the constants i/jj can be set 
to zero using an isometry of EI^"''^' ^° that we have left N{k + 1) parameters. 
These correspond to Nk charges and angular momenta. The N distances 
and the N masses are fixed by the choice of the intervals Ij. 

There is a partial converse to Theorem |l[ Let ip = {u,v,x,''P) '■ — > 
be a solution of the Reduced Problem. Let {p, (j), z) be cylindrical 
coordinates on M'^. li uj = 2{dv + x ' dip + ip ■ dx), then it follows from the 
harmonic map equations that the two-form e^"ig * u; is closed, and hence 
there is a function w such that dw = e^^i^ * lo. Here * is the Hodge star 
operator of the Euclidean metric on M^. Similarly, one checks that there is 
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a function A such that 
dX = du + p 



dp 



+ 2p 



UpU^ + ^e'^^'^p (^z + e^"(Xp ■Xz + tpp-'iPz 



dz. 



Also, there is an M'^-valued function 6 such that d9 = e^"i^ * dtf^ — w dx- 
As before, the forms dX and dO are closed thanks to the harmonic map 
equations. We now define the metric g on M' = M x (M^ \ A) by: 

and the M^-valued one-form A by: 

A = -{xd(t) + edt). 

We have: 

Theorem 2. Lei = {u,v,x,ip)- M'^ \ S ^ Hc^"*^ a solution of the 
Reduced Problem, and define M, g, and A as above. Then {A4,g,A) is a 
solution of the Einstein/ Ahelian- Yang- Mills Equations. 

The intervals Ij C A correspond to an event horizon in {M,g). However, 
as we mentioned in the introduction, unless = 1, it is expected that the 
metric g cannot be extended across all of S. The obstruction is a conical 
singularity on some bounded components of S. This singularity can be 
related to the force between the black holes, see [ pWj IWel - Nevertheless, 
(M, g) should be asymptotically flat, since a conical singularity cannot occur 
on the two unbounded components of S. To establish asymptotic flatness, 
the first step is a 'regularity' result for across S as in |Wel, LTi2]. This 
will be addressed in a future paper. 

3. Results and Proofs 
In this section, we prove the following theorem which is a generalization 



of the Main Theorem in |We4]. 



Main Theorem . Let H he a closed smooth suhmanifold ofW^ of codimen- 
sion at least two, and let denote its connected components. For each 
1 ^ j ^ N , let ifj'. \ S — 5- H 6e a T,-singular harmonic map into some 
geodesic 7j of the space H. Suppose there is a map if: M" \ S ^ IH which 
coincides in a neighborhood of each T,j with ipj, and for which there is a 
constant c > such that 

|r(^)| ^c(l+r2)-3/2. 



(5) 

Then, there is a unique harmonic map ip: 
to if near each Sj, and at infinity. 



which is asymptotic 
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In view of Example ^, we obtain as an immediate corollary the existence 
and uniqueness of solutions to the Reduced Problem. 

Corollary 1. The Reduced Problem for the stationary and axially symmet- 
ric Einstein/ Ahelian- Yang-Mills Equations has a unique solution. 



Proof of the Main Theorem. The existence proof is divided, as in | We4], into 
three steps: 

Step 1: for each sufficiently large ball Br C M", there exists a map 
'-PR ■ Br \ S ^ H such that ipR is asymptotic to ip near each T,j , and 
such that ipR = If on OBr. 

Step 2: there is a constant C independent of R such that dist((^/j, (p) ^ 
C in Bn\^. 

Step 3: there is a sequence -R, — > oo for which (pR. converge uniformly 
on compact subsets of M"" \ S to a harmonic map ip : M" \ S — > H which 
is asymptotic to (p near each T,j, and at infinity. 



Step 1. This follows from Proposition 1 in |We4]. Although, we required 
there that p was harmonic, this was never used in the proof of Proposition 1. 
We also required that dist((^(x), 7j) — > as x — Sj, but that is clearly 
satisfied here since (p agrees with pj near T,j. In fact, the existence of (pn is 
proved using a variational approach, and thus one also obtains in the course 
of the proof an energy estimate which is used in Step 3. 

We may assume, without loss of generality, that all the geodesies 7^ have a 
common initial point p = ^j(—oo) € dM. We use on H the coordinate system 
(u, d): EI ^ M X M"*"^ mentioned in Section |2|, where u is the Busemann 
function associated with the point p: 

u{q) = lim (dist {q, 71 (t)) + t) , 



see I We3| , Lemma 6]. The metric then takes the form (||). Write p = {u,v) 



and Pj = {uo,Vj), where uq is a harmonic function with uo{x) ^ 00 as 
— 5- S. Here, we have assumed that all the maps pj have the same harmonic 

function uq as their u coordinates, but that can easily be arranged. 

Let Hi{Br) be the Sobolev space of functions u on Br such that u, Vu G 

L'^{Br), and denote the closure of C^{Bji) in that space by Hifi{BR). Let 

Hf{BR,W"~^) be the weighted Sobolev space of R"*~"'^-valued functions v 

on Br such that 



{|t;p + Q^(V^)} < 00, 



'Br 

and let the closure of C^{Br \ S,R'"-i) in that space be HIq{Br,R'^-^] 
Now, define the space 'K{Br) of maps p = {u, v) on Br satisfying 

[dist{p,p)GL^{BR), 
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and let J{c(-B_r) be the space of maps (p G 'K{B^) for which dist(yj, f) ^ C. 
Define the renormahzed energy of = {u,v) € 'Kc{Bfi) by: 

Fr{v) = I I |V(^x - no)|' + Q^{Vv)\ . 

The solution is found by minimizing -Fr((^) over ip E 'K{Bji). In fact, an 
almost standard direct variational argument shows that a minimizer of 
exists in "KciBji) for any C > 0. The main difficulty is in proving that 
for some C = Cr large enough, depending on R, the infimum of over 
'^c{Br) is the same as over '^{{Br). 

Proposition 1. Let R > be large enough so that each bounded component 
of S is contained in the hall Br. Then there is a unique harmonic map 
$R - Br \ S ^ H such that (pR is asymptotic to sp near each Sj, and such 
that (fR = If on OBr. Furthermore, (pR € 'K{Br). 



Step 2. The main difference with | We4 | lies in this step. We must now 



show that the solutions <pR obtained in Step 1 belong to the space "He {Br) 
with the constant C independent of R. For this purpose, we will need two 
simple lemmas. 

Lemma 1. Let Lpi,{p2- ^ M. be smooth maps, and let p = dist((/9i, 992)- 
Then, we have 



(6) A\/l + p2^- lr((^i)| + |r(<^2)| 



Proof. Let Gp{q) denote the gradient of the function q 1-^ dist(p, g) on H 
evaluated at q ^ p, and let H denote the Hessian of the function (p, q) 1— > 
dist(p, q) on EI X H as a quadratic form on T{M x H) away from the diagonal 
{{p,p) G H X H}. The proof of the lemma relies on the following identity: 

Ap = (G^i(v32),t(v32))h + {G^^{pi),t{pi))m + H{dpi + dp2,dpi + dp2), 

which holds wherever p ^ 0. Equation (P) now follows from the fact that 
|Gp(g)| = 1 for all q ^ p, and the fact that H ^ 0, due to the negative 
curvature of H, see [^Y| , p. 368]. □ 
The next lemma is a simple generalization of the maximum principle for 
subharmonic functions. 

Lemma 2. Suppose a e C°°{Br \ S) is bounded and satisfies: 

Ao- ^ -c(l + r^)-^/^ on Br\^ 
cr = 0, on OBr. 

Then, there holds: 

(7) a < -, on Br. 

n — 2 
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pf Clearly, it suffices to prove the lemma with c = 1. Let v be the solution 
of 

Ai. = (l + r2)-3/2 
— 5- 0, as x — > oo. 

Note that such a solution exists, is radially symmetric, negative, bounded, 
and satisfies 

(8) M. = .(„) = __i_. 

Now, we have A((t + i^) ^ on Br \ S, and o" + is bounded. Thus, we 



can show, using cut-off functions near S, as in |We4, Lemma 7], that a + u 
is weakly subharmonic on all of Br, with finite Dirichlet integral over any 
ball Bri CC Br. Thus, using the maximum principle, we obtain: 

(9) a + v ^ v{R) ^ 0. 

In view of (^), this implies (0). □ 
Step 2 now follows easily, for if p = dist(^K,^), and a = ^/\ + — 1 
then (T is bounded, and cr = on BBr. Using (^), it follows from Lemma |^ 
that a satisfies: 

Aa^ -c(l + r2)V'- 

Thus, by Lemma |2|, there is a constant C, independent of R, such that 
p ^ C. In Step 3, we will also use the more precise estimate (M). 



Step 3. Again, the proof of this step is very similar to the one in [We4], 
and we only give a sketch. By Steps 1 and 2, we have for each R large 
enough, a harmonic map (pR = (ur,vr) £ "KciBR). Now fix Rq, and for 
R > Rq consider the maps ipR restricted to Br^. It is not difficult to check. 



as in [We4], that FR^^ipR) is uniformly bounded, i.e., bounded independently 
of R. In fact, we have A{ur — uq) ^ 0, and 

\ur — uo\ ^ dist{ipR, (f) + \u — uo\ ^ C + sup \u — uq \ , 
on Br^^ \ T,. Thus, ur — uq are uniformly bounded subharmonic functions 



on S/jg \ S, and the argument used in the proof of [We4, Lemma 7], shows 
that the Dirichlet integrals of ur — uq over Br^ are uniformly bounded. The 
bound on the second term in FR^^ipR) now follows from the inequality: 

2Q^^{Vvr) ^ A{uR-uo). 

Therefore, we can find a sequence -Rj — > oo such that ipR- converges pointwise 
a.e. in Br. By a standard diagonal argument, we can assume that the same 
sequence works for all Rq. It is not difficult to see that the pointwise limit 
(/? is a harmonic map on M" \ S. Indeed if CC M" \ S, then for i large 
enough, ipR. \q is a family of smooth harmonic maps with uniformly bounded 
energy which maps into a fixed compact set of H. Standard harmonic map 
theory now implies uniform bounds in C^'"(r2), hence some subsequence 
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converges uniformly, by necessity to together with two derivatives. Thus, 
(/? is harmonic. Clearly, dist((/?, (p) ^ C hence if is asymptotic to near each 
T,j. It remains to see that 93 is asymptotic to (p at infinity. This follows 
from (0) applied to 



(TR = + dist(^R,(/?)2 - 1. 

In the limit we obtain 

a = \/l + dist((/9, — 1 ^ — c I/, 

which implies that dist((/3, ^) ^ as x ^ 00 in R" \ S. 

The proof of the uniqueness statement is unchanged from | We4 |. □ 
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